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Abstract 

-- ■ We propose a modification of tlie Faddeev-Popov procedure to construct a path 

Qv^ , integral representation for the transition amplitude and the partition function for 

0> \ gauge theories whose orbit space has a non-Euclidean geometry. Our approach is 

based on the Kato- Trotter product formula modified appropriately to incorporate the 

gauge invariance condition, and thereby equivalence to the Dirac operator formalism 
.<^ \ is guaranteed by construction. The modified path integral provides a solution to the 

^^ ' Gribov obstruction as well as to the operator ordering problem when the orbit space 

^vq . has curvature. A few explicit examples are given to illustrate new features of the 

formalism developed. The method is applied to the Kogut-Susskind lattice gauge 
%^J I theory to develop a nonperturbative functional integral for a quantum Yang-Mills 

[-., ' theory. Feynman's conjecture about a relation between the mass gap and the orbit 

space geometry in gluodynamics is discussed in the framework of the modified path 

integral. 
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0\ , 1. Motivations. In what follows we consider only gauge theories of a special (Yang- 

Mills) type where gauge transformations are linear transformations in the total configura- 
tion space. Let A" be a configuration space which is assumed to be a Euclidean space M^ 
p ■ unless specified otherwise. It is a representation space of a gauge group Q so the action of 

1^ ■ ^ on A" is given by a linear transformation x —>■ Q{uj)x. A formal sum over paths for the 

action invariant under gauge transformations would diverge. To regularize it, Faddeev and 
Popov have proposed to insert the identity ^ 
>< 

.c?; i = AM^) I dfi{n)6{x{^x)) , (1) 



where xi^) = is a gauge condition and App is the corresponding Faddeev-Popov determi- 
nant, into the divergent path integral measure so that after a change of variables x -^ Qx 
the (divergent) volume of the gauge group can be factorized out. The determinant App is 
calculated by doing the group averaging integral (|l]). 

This formalism has provided a solution to two important problems in the perturbative 
path integral quantization of gauge theories: The unitarity of the perturbative S-matrix 
and the construction of a local effective gauge fixed action. Later on it has been observed 
that the procedure becomes ill-defined beyond perturbation theory [Q]. A geometrical rea- 
son for the Gribov obstruction is that there is no gauge condition which would provide a 
global parameterization of the gauge orbit space A4 = X/Q by a set of affine coordinates 
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without singularities [Q]. This situation is often rendered concrete in the vanishing or even 
sign changing of the Faddeev-Popov determinant. An additional problem expected in the 
nonperturbative regime is the operator ordering ambiguity [Q . When the kinetic energy x"^ 
is projected on the gauge orbit space A4, one gets {h,gphh), where h are local coordinates 
on Ai and Qph is an induced metric on Ai. The metric is not flat |^. In order to establish 
a correspondence between the Dirac operator formalism and the gauge-fixed path integral, 
the latter must be appropriately defined to take into account the operator ordering in the 
physical kinetic energy operator. It is known that the formal path integral for motion on 
a curved manifold is ambiguous 0. Up to the second order of perturbative Yang-Mills 
theory, operator ordering corrections have no effect on the physical S-matrix 0]. Their role 
in a nonperturbative description is unknown. A standard argument used in dimensional 
regularization, ihS^{0) = 0, (a commutator of canonically conjugated field variables at the 
same spacetime point, which typically emerges in the operator ordering terms) does not 
seem to be applicable because the validity of dimensional regularization beyond perturba- 
tion theory is not justified [0]. Finally, the space A4 may have a nontrivial topology [0]. 
Due to the locality of the Faddeev-Popov construction, global properties of A1 are lost in 
the path integral. On the other hand the topology of the orbit space may be important for 
the spectrum of physical excitations. Feynman conjectured that the mass gap in gluody- 
namics in (2+1) spacetime might be caused by a compactification of certain directions in 
the configuration space upon an identification of gauge equivalent configurations |^ . 

In this letter we propose a modification of the procedure (|l|) which gives a solution to 
the above problems and leads to a well-defined path integral on the orbit space M. Its 
validity is demonstrated with a few examples. 

2. A modified Kato-Trotter product formula. Let H and o"a be self-adjoint 
Hamiltonian and operators of constraints, respectively. We also assume [H, o"a] = which 
holds for a sufficiently large class of physically interesting gauge models. Let the gauge 
group be compact. Consider a projection operator 

V= f rf/i(cu)e*""'^'" , f dfiiu) = 1 , P = pt = p2 _ (2) 



The operator (H) projects the total Hilbert space onto the Dirac subspace of gauge invariant 
states. If the gauge group generated by aa is not compact, we take a sequence csVs of 
rescaled projection operators Vs, where Vs projects on the subspace X^'^a — '^^5 ^^^ then 
consider the limit 6 ^ in a weak sense. The coefficients cs can always be chosen so 
that the limit exists and defines a Hilbert space isomorphic to the Dirac one 0. Let the 
Hamiltonian have the form H = Hq + V, where V = V{x) is an everywhere regular potential 
and Hq = p'^/2 is a kinetic energy operator. The operators Ho and V are assumed to be 
gauge invariant, i.e., they commute with the constraint operators^. Consider the evolution 
operator Ut = exp{—iTH) in the total Hilbert space. The physical evolution operator is 

■^This assumption is justified for gauge theories of the Yang-Mills type. More general cases will be 
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then f/f* = UtV. Making use of the Kato-Trotter product formula |TT]| for the evolution 
operator Ut we get the following product formula for the gauge invariant evolution operator 



where e = T/n and n ^ oo. In our previous work [|T^] we have shown that in the continuum 



limit the averaging gauge variables u play the role of the Lagrange multipliers for the 
constraints in the classical action. By making use of the classical theory of Kolmogorov, 
we gave an explicit construction of a countably additive probabihty measure for the gauge 
functions uj{t). One natural choice was such that any set of values of uj{t) at any set of 
distinct times is equally likely. Now we show that this concept leads to a modification of 
the Faddeev-Popov procedure. 

First we observe that by construction 

U^{x, nx') = U^{nx, x') = U^{x, x') (4) 

that is, the projected amplitude determines a genuine transition amplitude on the orbit 
space A^. To obtain its relation to the conventional gauge fixed amplitude, we introduce a 
set of local coordinates h on Ai via a change of variables 

X = Q{uj)h , (5) 

where the configurations x = h satisfy a gauge condition x(x) = 0, i.e., h E H C X and 
without loss of generality we may assume Ti ~ M^'^ (unless specified otherwise) with M 
being the number of physical degrees of freedom. The gauge invariant amplitude (|^) can 
be reduced on the gauge fixing surface H: U^{x, x') = U^{h, h'). This leads to a necessary 
modification of the procedure (P: Instead of averaging the gauge condition over the gauge 
group, we propose to average an infinitesimal evolution operator kernel for a free motion 
as follows from the modified Kato-Trotter product formula (|^) 



Since e — *> 0, the averaging integral can be done to proper order in e explicitly by the 
stationary phase approximation. The kernel U^{h,h') is obtained by adding —ieV{h) to 
the exponential in (||) in accordance with (j^) because V{x) = V{h) thanks to the gauge 
invariance of V. 

Suppose there exist transformations of the new variables a; and h in (j^, u -^ u + Ug 
and h -^ hg, such that x{uj, h) = x{uj + ujs, h^), where all configurations Xg = hs also satisfy 
the gauge condition x{h) = x{hs) = 0. From Q{uj + ujs)hs = Q{uj)h follows that hg = Qgh, 
i.e., hs is a Gribov copy oi x = h on the gauge fixing surface x{^) =0. If S is a collection 
of all residual gauge transformations, the modular domain K is the quotient 7i/§. The 
relation (j^) determines a one-to-one correspondence ii h & K and u parameterizes a group 
manifold. If h, h' G K, then the following relation holds 

[/^(/i, h';) = U^{hs, h') = U^{h, h') . (7) 



It determines an S — invariant continuation of the transition amplitude to the covering 
space Ti of the modular domain K. The Faddeev-Popov determinant specifies the volume 
element on M. |^, so for any integral involving only functions on M., f{x) = f{Qx), one has 
J dxf{x) = Vg fj^dhApp{h)f{h), where Vq is the volume of the gauge group. Therefore 
the folding of two projected infinitesimal evolution operators reads 

UFeih, h') = [ dhAFp{h,)U^{h, h)U^ih,, h') . (8) 

JK 

The restriction of the integration domain in (|]) does not lead to a formal restriction of 
the integration domain in the path integral resulting from (^ in the continuum limit as 
one might naively expect from (^. An important point is that, ii K ^ 7i (the Gribov 
problem is present), then U^^{h, h') has no standard form of an infinitesimal free transition 
amplitude on a manifold proposed in as we now proceed to demonstrate. 

The stationary phase approximation can be applied before the reduction of UP{x,x') 
on a gauge fixing surface. A deviation from the conventional gauge-fixing procedure results 
from the fact that there may be more than just one stationary point. We can always shift the 
origin of the averaging variable u so that one of the stationary points is at the origin uj = 0. 
Let Sa be operators generating gauge transformations in X. Decomposing the distance (x — 
Q{uj)x'y in the vicinity of the stationary point, we find (x — x', 6ax') = 0. So in the formal 
continuum limit we get the condition o"a(x,x) = {x,dax) = induced by the averaging 
procedure. This is nothing but the Gauss law enforcement for trajectories contributing 
to Up{x,x'). Suppose there exists a gauge condition Xa{x) = 0, which involves no time 
derivatives, such that a generic configuration x = h satisfying it also fulfills identically the 
Gauss law, i.e., {h, Sah) = 0. We will call it a natural gauge. In this case all other stationary 
points in the integral (|^) are Uc = ujg where Q{uJs) = i^^ G §. Therefore we get a sum over 
the stationary points in the averaging integral (^) if the Gribov problem is present. Still, in 
the continuum limit we have to control all terms of order e. This means that we need not 
only the leading term in the stationary phase approximation of (^ but also the next two 
corrections to it. Therefore the group element Q{u!) should be decomposed up to order u^ 
because cj^/e ~ e as one is easily convinced by rescaling the integration variable u -^ \feijj. 
The measure should also be decomposed up to the necessary order to control the e-terms. 
The latter would yield quantum corrections to the classical potential associated with the 
operator ordering in the kinetic energy operator on the orbit space. We stress that the 
averaging procedure gives a unique ordering so that the integral is invariant under general 
coordinate transformations on A^, i.e., does not depend on the choice of x (cf. (|14D). Thus, 



Uf'ih.U) = (27rze)-^^/2^Z}-i/2^/i,/i;)p('^-'^^)'/2^-*^^^('^''^^) + 0(e^)] , (9) 

s 

= Y,D-'^\h,K)U,{h,K) , (10) 

where D(h, h') is the conventional determinant arising in the stationary phase approxima- 
tion, and by Vq we denote a contribution of all relevant corrections to the leading order. 
The amplitude U['{h,h') is obtained by adding —ieV{h) to the exponential in (|). Note 
that V{h) = V{hs) thanks to the gauge invariance of the potential. 



In general, the equations aa{x,x) = are not integrable therefore the functions Xa 
discussed above do not exist. In this case we consider two possibilities. Let fic{h, h') be the 
group element at a stationary point in (^). Decomposing the distance in the vicinity of the 
stationary point we get [h—VLji'-, Sa^ch') = 0. Let x be such that the latter condition is also 
satisfied if h' is replaced by h'^. In this case the sum over the stationary points is again given 
by (^) where in the first term of the exponential h'^ is replaced by Qc{h, h'^)h'g. In the most 
general case, the sum over stationary points does not coincide with the sum over Gribov 
copies. However for sufficiently small e, the averaged short-time transition amplitude can 
always be represented in the form ([l0|). Note that as e goes to zero, f/°(x, x') -^ 6'^{x — x'). 
If we change the variables in the delta function according to (||) and then do the gauge group 
average, we get, in general, a sum of delta functions on the gauge fixing surface. From the 
symmetry of the change of variables (|^) it follows that the support of the averaged delta 
function is formed by points h = h'^, h' & K and h E Ti, which means that the sum over 
S emerges again for sufficiently small e. We will explain below how to calculate relevant 
corrections to the leading order of the limit e — *> 0. 

Let us replace U^{h,hi) in (H) by the sum (p!OD and make use of (0) applied to the 
second kernel in (^: U^{hi,h') = U^{his,h'). Since the measure on the orbit space does 
not depend on a particular choice of the modular domain, dhs^Fp{hs) = dhApplh), we 
can extend the integration to the entire covering space by removing the sum over S 

U^eiK h') = j dh\AFp{h)\D~'/\h, h)U,ih, h)U^ih, h') . (11) 

The absolute value bars accounts for a possible sign change of the Faddeev-Popov deter- 
minant on the gauge fixing surface. The procedure can be repeated from left to right in 
the folding (^ , thus removing the restriction of the integration domain and the sum over 
copies in all intermediate times t G (0,T). The sum over S for the initial configuration h' 
remains in the integral. 

Now we can formally take a continuum limit with the result 

h{T)=h 
U^(h, h') = J2 [AFp{h)AFpiK)r'^^ I Vhy^d^t^,e'-fodii(>^'9r,Hh)/2-V,{h)-V{h)] ^ (^2) 



where gph is the induced metric on Ai. The local density Ylt A/det gph should be understood 
as the result of the integration over the momenta in the corresponding time-sliced phase- 
space path integral where the kinetic energy is {ph,gphPh)/'^ with ph being a canonical 
momentum for h. To obtain (0), one usually sets h' = h — A for h = h{t) and h' = 
h{t — e) in each intermediate moment of time t and makes a decomposition into the power 
series over A. According to the relation, obtained in ||^, between the volume of a gauge 
orbit through x = h, the induced metric Qph, and the Faddeev-Popov determinant we get 
D{h' + A, h') = A^pp{h')/ det Qphih') + 0(A). This relation explains the cancellation of the 
absolute value of the Faddeev-Popov determinant in the folding (^ made in accordance 
with the rule ([TT|). The term A^/e in the exponential (^ gives rise to the kinetic energy 
(A,^,,A)/2e + 0(A3). 



The most difficult part to calculate are the operator ordering corrections Vq{h) in the 
continuum limit. Here we remark that D{h,h') has to be decomposed up to order A^, 
while the exponential in (||) up to order A"^ because the measure has support on paths 
for which A^ ~ e (i.e., A^/e ~ e). There is a technique, called the equivalence rules for 
Lagrangian path integrals on manifolds, which allows one to convert terms A^" into terms 
e" and thereby to calculate V„ P, ^2 



p{iivj2fe) 

where the sum is extended over all permutations of the indices j to make the right-hand side 
of ( |13|) symmetric under permutations of the j's. Following ([13| ) one derives the Schrodinger 
equation for the physical amlitude (0). The corresponding Hamiltonian operator on the 
orbit space has the form 

Hph = -l/2A]^],a, {[g;^y'AFpdk) + V{h) , (14) 

where dj = d/dhj. Observe that the kinetic energy in (0) does not coincide with the 
Laplace-Beltrami operator on Ai because [det gph]^^'^ 7^ App as shown in f^. The operator 
( [T^ ) is hermitian with respect to the scalar product where the volume element is dhApp. It 
is also invariant under general coordinate transformations on Ai, i.e., its spectrum does not 
depend on the choice of local coordinates on M. and, therefore, is gauge invariant. We skip 
technical details of the derivation of ([T^) because they are standard in the path integral 
formalism and instead turn to examples to illustrate the main features of the new projection 
formalism. But before we do so, let us comment that the effects on the energy spectrum 
caused by the modification of the path integral can be found from the pole structure of the 
trace of the resolvent 

tr^(r) = ti (r-iHY^ = f dTe-^^tiU^, (15) 

trf/^ = I dhApp{h)U^{h,h) . (16) 

Jk 

3. Examples. Let us take the simplest example ^ where x G M'^ and the gauge transfor- 
mations are rotations of x about the origin. The orbit space is a space of concentric spheres. 
The natural gauge is x* = 6^^h. This gauge is not complete because there are rotations that 
change the sign oi h: h —>■ ±h. These transformations form the group § = Z2. The mod- 
ular domain is an open half-axis K = (0, 00). The measure on the orbit space is obtained 
by introducing spherical coordinates on M^ so that h = \x\, while the spherical angles pa- 
rameterize the gauge group manifold SO{3)/SO{2). It is noteworthy to observe that in 
the projection formalism there is no principal difference between reducible and irreducible 
gauge theories. Here we obviously have a reducible case because x has a stationary group 
SO {2). The Faddeev- Popov determinant is defined with respect to a subset of independent 
constraints. In this case, the constraints are three components of the angular momentum. 
So we can take the angular momenta about the axes i = 2, 3 as independent constraints and 
find that App = h"^ in the full accordance with the measure resulting from the spherical 
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coordinates. The angular measure must be normalized by Ait to fulfill the condition in (j^). 
So Eq. (H) assumes the form 

U^^{h, h') = (27rie)-=^/2 1/2 F dujsmiue'^^'+''"^/^'-'^''' ""'^Z' . (17) 

Jo 

Doing the averaging integral exactly and substituting the result into @ we obtain 

U^{h, h') = [hh']-^ {Urih, h') - Urih, -h')} , (18) 

where the amplitude UT{h,h') is given by a standard path integral for a one- dimensional 
system with the gauge-fixed action Sgf[h] = J^ dt[h'^/2 — V{h)]. There is no operator 
ordering correction. 

The stationary phase approximation leads to the same result if one takes into account 
all the stationary points, i.e., Ug = and Ug = vr. The second stationary point is associated 
with the only nontrivial residual gauge transformation h' — > —h'. Here one also gets Vq = 
(cf. (^). The same path integral can be derived from the Schrodinger equation projected 
on the orbit space [|l^. It is not hard to calculate the partition function (0) for a harmonic 



oscillator V = x^/2 = /i^/2 and find the resolvent (p!5|). The poles determine the spectrum 



En = 2n + 3/2. The distance between the energy levels is 2, while the oscillator frequency 
in the Hamiltonian is 1. The effect of the doubling would be lost, had we ignored the 
redundant discrete gauge symmetry and calculated the partition function by a standard 
path integral for the gauge- fixed action Sgf. 

A simple model to illustrate the effect of the residual gauge symmetry in systems with 
several physical degrees of freedom is a Yang-Mills theory in (1+0) spacetime [0. A 
total configuration space is a Lie algebra X, and the gauge group acts in the adjoint 
representation x — > QxQ~^. The quadratic form in the exponential (|^) is taken with 
respect to the Killing form on X: {x,x) = trx^, where xy = [x,y] for any x,y & X and 
the commutator stands for the Lie bracket in X. If a matrix representation is assumed, for 
compact groups it can always be normalized to the ordinary matrix trace. The Gauss law 
is cr = [x, x] = 0, therefore the natural gauge condition to parameterize A4 is x = h E H 
where if is a Cartan subalgebra (a maximal commutative subalgebra of X) . The system 
has r = dim if physical degrees of freedom. The invariant measure d^{uj) is proportional 
to (ia;det[(expa; — l)u)~^]. In the exponential in (H) we get {h — exp(a;)/i')^. An adjoint 
transformation of x to if is not unique and determined modulo the Weyl group W acting in 
H. The Weyl transformations are inequivalent compositions of reflections in the hyperplanes 
orthogonal to simple roots in H |T^. So S = W and M. ~ H/W = K'^, where K~^ is the 



Weyl chamber in H. The stationary points of the integral (|g) are cUc = uJw where the 
adjoint action of the group element expcuvy on h transforms the latter by a Weyl reflection. 
Introducing the Cartan- Weyl basis [^ one can find the Jacobian of the change of variables 
(|^) (the Faddeev-Popov determinant) App{h) = deth = K^iji) where n{h) = na>o('^'^)' 
the product is extended over all positive roots in H. The operator h acts in the orthogonal 
complement of the Cartan subalgebra. Note that the Cartan group is the stationary group 
oi h E H, so uj should take its values in the orthogonal supplement of the Cartan subalgebra. 
The factor D{h, h') is given by the determinant of the operator hh'. Doing the stationary 
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phase approximation in the averaging integral, and then taking the folding (j^) we find 

U^{h, h') = Y,Hh)K{h'^)]-^UT{h, h'^) . (19) 

w 

Here the amplitude Uxih, h') is given by the conventional path integral (no restriction of the 
integration domain to the Weyl chamber) with the action j dt{h? /2 — V{h)). There is no 
operator ordering correction, and the physical metric is Euclidean. The same integral follows 
from a direct solution of the Schrodinger equation [|I^ . The Faddeev-Popov determinant 
vanishes at the hyperplanes (/i, a) = 0. The boundary of the modular domain (the Weyl 
chamber) is formed by these hyperplanes when a runs over the simple roots. The amplitude 
([l9| ) remains finite at the singular points (the Gribov horizon) because K,{hw) = Pwi^{h) 
where pw = ±1 is the parity of the corresponding Weyl transformation (cf. also (|TB|): It 
remains finite at h = or h' = 0)). Contributions of paths hitting the horizon have a 
dramatic effect on the spectrum of the isotropic oscillator V = x'^/2 = h'^/2. Doing the 
path integral for Uxih^h') in (|1^) and calculating the poles of the resolvent (0) we find 
llTSl En = viTii + ■ ■ ■ + Vr^r + N/2, whcrc Tij are non-negative integers, z/j are orders of 
invariant irreducible symmetric tensors of the Lie algebra, and r and A^ are, respectively, 
the rank and dimension of the group (e.g., for SU(r+l), Ui = i + 1, i = 1,2, ...,r). Thus, the 
physical frequencies of the isotropic oscillator appear to he Ui y^ 1. The partition function 
calculated by a standard path integral for the gauge fixed action of the model would lead 
to the spectrum of the isotropic oscillator (i.e., z/j — > 1). 

To illustrate the effects of curvature of the orbit space, we consider a simple gauge 
matrix model. Let x be a real 2x2 matrix subject to the gauge transformations x — > Q{u})x 
where Q G S0{2). An invariant scalar product reads {x,x') = trx^x' with x^ being a 
transposed matrix x. The total configuration space is M^. Let e be a generator of S0(2), 
i.e., £ij = —£ji, and £12 = 1. Then Q{uj) = exp^ue). The Gauss law enforced by the 
projection is a = trx^ex = 0. It is not integrable. We parameterize A^ by triangular 
matrices /i2i = (the gauge 0:21 = 0). The residual gauge transformations form the group 
S = Z2: h -^ ±h. The modular domain is a positive half-space hn > 0. Calculating 
the Jacobian of the change of variables (|^) for this model we find the Faddeev-Popov 
determinant App{h) = hn. The plane hu = is the Gribov horizon. The averaging 
measure in (|^) reads {27r)~'^duj, where uj G [0,27r). The quadratic form in the exponential 
in (H) is 

[h-e'^^h'Y = {h,h) + {h',h')-2{h,h')cosuj~2{h,sh')smuj . (20) 

A distinguished feature of this model from those considered above is that the stationary 
point is a function of h and h'. Taking the derivative of (|20|) with respect to u and setting 
it to zero, we find 

1 (h,eh') , , , 

ujc = tan-i ,; ' , uj', = Uc + n . 21 

[h, h') 

The second stationary point uj^ is associated with the Gribov transformation h -^ —h. A 
geometrical meaning of the transformation h' — > exY>{ujc£)h' is transparent. The distance 
[{h — h'YY/"^ between two points on the gauge fixing plane is greater than the minimal 
distance between the two gauge orbits through x = h and x' = h'. By shifting x' along the 



gauge orbit to x'^ = exp{ujc£)h' a minimum of the distance between the orbits is achieved. 
In such a way the metric on the orbit space emerges in the projection formahsm. Its explicit 
form can be found as has been explained in section 2. We substitute u = ujc{h, h') into 
([20|), set h' = h — A and decompose (^) in a power series over A. The quadratic term (the 
leading term) determines the metric. We get {A,gph{h)A) = {A,A) + {A,£h){£h,A)/{h,h). 
The metric is not flat. The scalar Riemann curvature oi A4 is R = 6/{h,h) = 6/{x,x). 
Note that it is gauge invariant as it should be because the curvature R does not depend on 
the choice of local coordinates on A^. In the stationary phase approximation the cosine and 



sine in (^^ should be decomposed up to fourth order in the vicinity of the stationary point. 
In this model quantum corrections do not vanish. The infinitesimal transition amplitude 
on Ai assumes the form 

U^{h,h') = (27r2e)-=^/2r^-i/2(^^^/^giS.(M')^^-i/2(^^_^/)gi5.(h,-/^')l ^ (22) 

S,{h,h') = [{h,h) + {h',h')-2D{h,h')]/{2e)~e[8D{h,h')]-^-eV{h) , (23) 

where D{h, h') = {h, h') cosUc + {h,eh') smujc- Observe that the middle term in (P3|) is a 



quantum correction to the classical potential in the short-time action (p3|) (it is proportional 
to h"^, if one restores the Planck constant). The two first terms in (^) determine [uniquelyl) 
the operator ordering in the kinetic energy operator on the orbit space in accordance with 

Finally, we consider the case of an infinite number of stationary points in the averaging 
integral (^. This would have an effect that the modular domain K may become compact, 
which, in turn, leads to a discrete spectrum regardless of the details of the potential and 
therefore gives a mechanism for the gap between the ground state and the first excited state. 
A suitable soluble model is the Yang-Mills theory in cylindrical space time M, ITR IIBI where 



space is compactified into a circle. A Yang-Mills theory without matter in two dimensions 
has no physical degrees of freedom unless the spacetime has a nontrivial topology. The latter 
is exactly our case. The total configuration space X is the space of all periodic connections 
A{x + 27c) = A{x) taking their values in the Lie algebra of a compact Lie group. The gauge 
group Q acts on X a.s A ^ ^A = QAQ~^ — iQdQ~^, and Q{x + 2tt) = Q{x) (i.e., the gauge 
group elements Q can be continuously deformed towards the group unit [^, 0). The 
exponential in (^) assumes the form J '" dxtT{A — ^A')^ for any two configurations A{x) and 
A'{x). It is the distance between A{x) and ^A'{x) introduced by Feynman [Q. A rigorous 
definition of the averaging integral over Q can be given within the Kogut-Susskind lattice 
gauge theory, to which we turn shortly. Note, for example, that N in (^ is infinite because 
the number of degrees of freedom is determined by the number of Fourier modes of A{x). 
Here we assume the existence of a normalized averaging measure on Q. The gauge 
group average enforces the Gauss law a = dA — i[A, A] = V{A)A = 0. The orbit space A4 
can be parameterized by constant connections A{x) = h taking their values in the Cartan 
subalgebra H [|14|, |16|. This is the natural gauge because a{h, h) = 0. The number of the 
physical degrees of freedom equals the rank of the gauge group. We set Q{x) = exp{iu{x)). 
The averaging function u takes its values in the Lie algebra and its constant (independent 
of x) Fourier mode does not have any component in the Cartan subalgebra because the 
connection A{x) = h is invariant under constant gauge transformations from the Cartan 
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subgroup. The Gaussian integral over such functions can formally be done in (j^) yielding 
D~^/'^{h^h') = det ~^/^[— V(/i)V(/i')] which is easy to compute in the Cartan-Weyl basis 



16| . In the lattice formulation described below this procedure would not involve any 
infinite factors in D{h, h') as in the formal continuum limit. Here we draw attention to the 
following feature. A projection of any connection A{x) into the subspace of constant Cartan 
connections h is not unique and determined modulo the afiine Weyl transformations []I^]. 
The affine Weyl group W^ is a semidirect product of the Weyl group W and the group of 



translation along the group unit lattice in the Cartan subalgebra |TP[. Thus, the number 
of stationary points in the averaging integral @ is infinite because the number of elements 
Wa is infinite: h ^ hg = hw + 2Q;na/(a, a), where hy/ is a Weyl transformation of h, a 
any positive root, n^ an integer. The modular domain K = H/Wa is called the Weyl cell 
ll9| . It is compact in H. This has a significant effect on the spectrum of the system. In the 
continuum limit the transition amplitude has the form (^) where the sum is extended over 
the affine Weyl group (over an infinite number of stationary points /i'^)|l^, the Faddeev- 



Popov determinant reads App{h) = K'^{h) and K,{h) = nQ>o sin(/i, a). There is no potential 
energy in 1 + 1 dimensional Yang-Mills theory {V{h) = in ([19D). Thus, we get an amplitude 
for a free particle symmetrized over the affine Weyl group. 

In the simplest SU(2) case, the system has only one physical degree of freedom. The 
Weyl transformations are reflections hw = ^h and the translations over the group unit 
lattice are shifts h ^ h + 2n. The amplitude (|19D appears similar to the one for a particle 
in an inflnite well or on a circle. The spectrum of such a system is discrete. In general, 
we can calculate the poles of the resolvent (^) for the 2D Yang-Mills theory using our 
path integral formalism and flnd [|l^ the spectrum E^ = EoK'jn, In) — {p, p)] where Eq is a 



constant, p is a half the sum of all positive roots, and 7„ ranges over a lattice in the Weyl 
chamber that labels all irreducible representation of the gauge group. 

Thus, an inflnite number of stationary points in the averaging integral (j^) is, in fact, 
evidence that the physical conflguration space (or its subspace) is compact. This is the 
mathematical structure of the path integral which is to be sought when studying the relation 
between the mass gap in gluodynamics and the compactness of the orbit space. Physically, 
the existence of an inflnite number of the stationary points separated by different distances 
(large and small) implies that physical eigenstates of the Hamiltonian should exhibit a 
periodicity in some directions in X, which may cause the mass gap as has been argued by 
Feynman [Q. 

4. Kogut-Susskind lattice gauge theory. To fulflU our program in realistic Yang- 
Mills theory, we need a regularization at short distances to give a rigorous meaning to the 
averaging integral (^ which is a functional integral. The simplest possibility is to set the 
system on a lattice. Since we still want to have a Hamiltonian formalism, the time should 
remain continuous. With such a choice of the regularization we arrive at the Kogut-Susskind 
lattice gauge theory [^ . 

Let y be a three-vector whose components are integer-valued in the lattice spacing 
unit a, i.e., y labels the lattice sites. With each link connecting the lattice sites y and 
y + ia, where i is the unit vector in the direction of the ith coordinate axis, we associate a 
group element Uy^i called the link variable. The collection of all values of all link variables 



10 



forms the configuration space A" of the system. The gauge group Q acts on X by the rule 
Uy^i -^ QyUy^iQyl^^. Thcrcforc the averaging measure is a product of the normahzed Haar 
measures at each site: d^ = Y[y dfiH{Qy). The Kogut-Susskind Hamiltonian is H = Hq + V 
where the potential V is nothing but the color magnetic energy on the lattice. The kinetic 
energy is the sum of kinetic energies of a rigid rotators associated with each link variable 

Ho = l/2j2^Tiuy,u-}f . (24) 

An important observation is that the coupling of the rotators occurs only through the 
potential V. Therefore the "free" amplitude Uj^{{u} , {u'}) , where {u} is a collection of 
configurations of rotators Uy^i, is factorized into a product of the amplitudes for each rotator 
Uj'{uy^i,u' j). An explicit form of the latter is well-known [^ (a transition amplitude of 
a particle moving on a group manifold). This replaces the free motion amplitude in a 
Euclidean space in the averaging integral (^). The averaging integral 

[/^ ({^}, W}) = / n [df^Hi^y)] n f^°(^.., ^yU'y,ny+a^) (25) 

y y,i 

can be calculated via both the stationary phase approximation and the decomposition over 
the characters of the irreducible representations of the amplitudes U^{u,u') proposed in 



||21|| . Thus, we have reduced the problem of constructing the transition amplitude on the 
orbit space to a standard mathematical problem in group theory. It is important that the 
effects of the Riemann structure (metric) and topology of the orbit space are completely 
determined by the averaging procedure proposed and do not depend on the details of the 
self-interaction V. Note that the average in ( pS] ) induces a nontrivial coupling between the 
rotators because six link variables attached to the same site are transformed with the same 
group element. So the factorization of the amplitude on M. is no longer possible. This is 
the result of the enforcement of the Gauss law in the physical amplitude. 

Thus, we get, in principle, a constructive and consistent procedure to study the effects 
of the non-Euclidean geometry of the orbit space on physical quantities like the spectrum 
of low energy excitations. In particular, in the strong coupling limit {g large) the potential 



V can be treated by perturbation theory |20]. To introduce an explicit parameterization 



of the orbit space, one can, for example, make use of Morse theory for the Kogut-Susskind 
lattice gauge theory as proposed in [^. In this case, the initial configuration {u'} in the 



averaging integral (|25D is taken from the corresponding modular domain. The projected 
amplitude has a natural gauge invariant continuation to the covering space of the modular 
domain. 

5. Conclusions. We have proposed a self-consistent path mtegra/ quantization of gauge 
theories which resolves the Gribov obstruction and provides control for the operator ordering 
problem, and it does not rely explicitly on the Schrodinger equation on the orbit space. This 
was our main goal. The continuum limit in the Kogut-Susskind lattice gauge theory might 
still be a difficult analytical problem. However, it is believed that the techniques we have 
developed might be useful to study the effects of the non- Euclidean geometry and topology 
of the orbit space in some physically reasonable approximations. We have also demonstrated 
its effectiveness with several examples. 
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